Magnetism and superconductivity in underscreened Kondo chains 
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We present a one dimensional model of electrons coupled to localized moments of spin S > 1 in 
which magnetism and superconductivity interplay in a nontrivial manner. This model has a non- 
Fermi liquid ground state of the chiral spin liquid type. A non-conventional odd-frequency pairing 
is shown to be the dominant instability of the system, together with antiferromagnetism of the local 
moments. We argue that this model captures the physics of the Kondo-Heisenberg spin S — 1 chain, 
in the limit of strong Kondo coupling. Finally, we discuss briefly the effect of interchain coupling. 



PACS: 71.27.+a 71.10.Pm 74.20.Mn 75.30.Mb 



The interplay of magnetism and superconductivity is a 
fundamental problem in condensed matter physics. Su- 
perconductivity is associated with the pairing of electron 
states related by time reversal. Magnetic states, in which 
time reversal symmetry is lost, should therefore strongly 
compete with superconductivity. Thus it came as a sur- 
prise when it was discovered in 1984 by Schlabitz et al. 
[EJ that magnetism and superconductivity actually coex- 
isted in the heavy fermion compound URu2Si2. Since 
then, other heavy fermion superconductors were shown 
to present magnetic moments in their superconducting 
phase H . All these compounds contain rare earth or ac- 
tinide ions with very localized 4/ or 5/ orbitals, strongly 
interacting with the conduction band. (This is in contrast 
to such cases as the Chevrel phases where magnetism and 
superconductivity coexist because the magnetic moments 
responsible for magnetism are only very weakly coupled 
with the electrons that form the condensate ||.) 

The physics of heavy fermion compounds is believed to 
be described by the Kondo Lattice model in which con- 
duction electrons interact with the local moments asso- 
ciated with the localized / electrons, and a large amount 
of theoretical work was carried out. However, although 
the single impurity Kondo problem is well understood j| , 
the theoretical analysis of the Kondo Lattice model has 
proven extremely difficult. This is because the Kondo 
effect (the quenching of the local moments) competes, in 
the lattice problem, with the Rudermann-Kittel-Kasuya- 
Yosida interaction which orders the local moments. This 
frustration is believed to be at the origin of the rich 
phase diagram of heavy fermions systems. There exist at 
present various theories of the superconducting ground 
state of the Kondo Lattice Model. Conventional scenar- 
ios involve the pairing of fermionic quasiparticles into 
either a spin triplet "p-wave" state or a spin singlet "d- 
wave" state |||. A less conventional pairing j|] involves 
the formation of a spin singlet, isotropic, odd in time su- 
perconducting order parameter, known as odd-frequency 
pairing |||. 

In one dimension, the theoretical situation is more fa- 



vorable j9j , since there are powerful methods to deal with 
strong interactions. Bosonization jlo| and Density Ma- 
trix Renormalization Group Jll]] have clarified much of 
the physics of the single channel 5 = 1/2 Kondo lattice. 
One finds a paramagnetic metallic phase for small Kondo 
coupling and away from half filling. However, conven- 
tional superconducting fluctuations are strongly reduced 
|l2| . Adding a direct Heisenberg interaction between the 
spins has been shown to cause an enhancement of odd- 
frequency pairing correlations. However, the formation 
of a spin gap precluded the observation of fluctuations 
towards magnetic ordering p3| . 

The aim of this paper is to discuss a one dimensional 
Kondo Lattice model with underscreened moments in 
which magnetism can be expected to dominate. We will 
show that in this one dimensional problem strong fluctua- 
tions to composite superconductivity coexist with strong 
fluctuations towards magnetic ordering. 

The generalized Kondo lattice Hamiltonian in one di- 
mension is: 

H = -^( c L Cl+1 > CT + 4+1,^) 

1,(7 

+ J2f0i-S i+1 ) (1) 

i i 

where Cj iCr annihilates an electron, Si is a spin-S* operator, 
and the function / describes the spin-spin interaction. In 
the following, we will consider the case Xk "C /, t. 

We shall be interested in magnetism and underscreen- 
ing coexisting with superconductivity. To that purpose 
choose S > 1. Several forms of magnetic interaction need 
to be considered in this case. The simplest possibility is 
to take a Heisenberg coupling /(£$ • Si+±) = XhSi • Si+i, 
but no interesting magnetic effects would ensue. We 
would get at Xk = cither a spin gap for integer S 
or an effective spin 1/2 chain for half-odd integer S [ pd[ . 
Turning on Xk <C t,Xn, the case of half-odd integer S 
would reduce to exact screening, and a metal with a spin 
gap will form |^). The case of integer S would lead 
to a completely trivial result, the local moments being 
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already completely screened by the formation of the Hal- 
dane gap preventing the Kondo effect and leaving the 
electrons essentially free. 

Magnetic effects, on the other hand, would result 
from interactions in the chain, /, which lead to low 
energy dynamics described by effective Hamiltonians 
richer than the SU(2)i Wess-Zumino-Novikov-Witten 
(WZNW) model which govern the half-odd integer 5 
Heisenberg M discussed above. We shall discuss sys- 
tems with fixed points described by SU(2)2s (5 > 1). 
An example is provided by the integrable spin-S chains 
@H?], where / = P s = £jLi«j-Pj> Pj is the spin- 
j projector and aj = 2i=i* _1 - Such SU2s(2) models 
arise naturally under some circumstances: it is known 
in particular that the Heisenberg spin-1 chain can be 
described as a SU{2) 2 WZNW model perturbed by a 
mass term of the order of the Haldane gap. The 
generalization of this result to arbitrary spin 5 chains, 
based on a perturbed SU(2) 2 s WZNW model, was ob- 
tained in Ref. (2l|, (and will be further discussed be- 
low). Thus, our study should describe the Underscreened 
Kondo-Heisenberg lattice when the Kondo coupling ex- 
ceeds the Haldane gap. For 5 = 1, this corresponds to 
OAXh < Xk < Xh- We shall study the model away from 
half filling and for all 5, revealing the appearance of a 
critical point describing a non Fermi-liquid where mag- 
netism and superconductivity interplay p2]| . 

We now proceed to determine the low energy behavior 
of the theory. Consider the Hamiltonian in the contin- 
uum limit. According to the standard prescriptions of 



non-abelian bosonization [|10],^3 24 1, the electrons are de- 
scribed by the following Hamiltonian: 



H e \ = H p + H a 



h p = J dx^- [(rfi p ) 2 + [d-Apf] 

Ha = — - — / dx(J R .J R + Jl-Jl) 



(2) 



(3) 



where the canonically conjugate fields H p , 4> p describe the 
charge excitations, and the non abelian SU{2)\ currents 
Jr,l describe the electron spin excitations. The local 
moments are described in the low energy regime by the 
Hamiltonian: 



H„ 



2ttv s 



2(1 + 5) 



dx(S R .S R + Sl-Sl) (4) 



where S R , L are SU(2) 2 s WZNW currents. The Kondo 
interaction Xk at incommensurate filling becomes: 



X K a / dx(J R + J L )-{S L + S R )(x) 



(5) 



and preserves spin-charge separation. Carrying out stan- 
dard RG calculations we find that this interaction is a 
combination of terms that are (marginally) relevant in 



the RG sense and purely marginal terms. The former 
drive us to a strong coupling fixed point and we need 
a non-perturbative way to determine it. To do so we 
may neglect the marginal couplings J R .S R and Jl-Sl as 
well as the velocity difference between electron and lo- 
cal moments spin excitations, and check, subsequently, 
their relevance at the fixed point. This leads us to the 
following spin Hamiltonian: 



H 



spin — H i + H2 

Hi= dx 



(6) 



Ho 



dx 



-; 2nv - - a t 

2(1 -f S) ^ — 7^- Jl-Jl + Xko,S r .Jl 

2(1 -|- S) ^ L '^ L ^ — t^-Jr-Jr + XkclSl-Jh 



A similar Hamiltonian was proposed in Refs. [|25|,[26| to 
describe the spin sector of a Hubbard chain coupled to N 
spin 1/2 chains (chain cylinder model). The N — 2 fixed 
point was analyzed in detail using an exact solution at a 
Toulouse point [p5| , and the correlation functions of the 
composite order parameters were obtained f26f| . 

The fixed point Hamiltonian can be determined by ar- 
guments of chiral stabilization j27j, since Hi and H2 in 
Eq. (0) are both chirally asymmetric (only their sum 
is chirally symmetric) p8|j . We find, that the electron 
spin degrees of freedom are described by the coset CFT 
- SU ^su( 2 U )2's >2S ~ 1 ' or ecm i va l en tly> the minimal model- 
M2S-1, of central charge c = 1 — (2S+i)(2S+2) > w h ue 
the local moments spin excitations are described by the 
SU(2) 2S -i WZNW model: 



H* = 



SU(2) 1 ^SU(2) 2 s- 



517(2) 



2S 



SU{2) 



2S—1- 



(7) 



This fixed point describes an interesting interplay of mag- 
netism and superconductivity. We shall discuss the struc- 
ture of the ground state, the thermodynamics and then 
the correlation functions of the model. 

The ground state is a coset singlet, formed between 



electrons and local moments. A fraction 



of 



(2S+l)(2S+2) 

the electron spin degrees of freedom is absorbed by the 
spins, leading to its complete screening. This manifests 
itself in the magnetic susceptibility and the specific heat: 



C 



3(25) 
25 + 2 
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(25+l)(25 + 2) 



where in the latter quantity we included also charge 
contributions. This leads to a Wilson ratio: Rw = 

( 2S+?)fs+ 9 i)(2g-i) going to zero as S ~> 00 ' 
We now calculated the long distance behavior of the 

physical correlation functions. To do so we have to ex- 
press the physical operators in terms of the operators 
around the fixed point, which we proceed to discuss. The 
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primary operators of SU (2) jv model are <J>^ (j < £ ) and 
carry spin-j. The coset primaries, <fp?„ , (0 < j < 1/2, 
< f < S - 1/2 and < j" < S), are related to the 
SU(2)n WZNW spin-j primaries via the decomposition: 



V L,1 V L,2S-1 — 



E (8) 

li— j' I <i" < j'+i' 



The conformal weight of the coset primary is such that 
the two sides of the equality have the same conformal 
weight, implying that j" has left conformal weight 
fili+H + 1^+11 _ ^g±il) . Similar identities hold for 

the right component. Next, note that at the fixed point 
products of left primary operators of SU(2)\ and right 
SU(2)2s WZNW models have to decompose into a sum of 
products of primary operators of the left minimal model 
and right SU (2) 2 s-i model. Moreover, the total spin has 
to be conserved. Therefore, one can write the decompo- 
sition: 



J2 <;" l| '«.^ : 

\j-j"\<j'<j+j" 



(9) 



This decomposition satisfies the requirements on the in- 
dices in (j??,, . It is formally similar to a Clebsch-Gordan 
decomposition. 

Thus, at the fixed point, the spin operator is given by: 



S, 



n ~ a(S' R + S' L ){na) + ^Tr^g'Xna) 



(10) 



where S R L and g' are respectively the currents and the 
SU{2) matrix field of the SU{2) 2S -i WZNW model, and 
(t\'n 2 is the field of the minimal model with scaling di- 
mension 2(25+1^(25+2) • ^ S a resu ^! the spin-spin corre- 
lation function behaves as: 



(S n .S m ) 



(_)»- 



(n — m) 2 (n — m) s 



(11) 



where 5 s = (25+1^5+2) ^ ^ There is a strong tendency 
to antiferromagnetism. 



The fermion operator ip^, is given by: 



^(e p +4, p )/V2:i/2.i/2, 



J L 



9r 



(12) 



,1/2,1/2 

J L 



is the antiholo- 



where 6 p (x) = tt f x U p (x')dx', 

1/21/2 

morphic component of the <p ' field of the minimal 
model, and g' R is the holomorphic component of the g' 
field of the SU(2) 2 s-i WZNW model. A similar expres- 
sion holds for ipR, in which L and R are exchanged and 
<j) p — > — p . The resulting fermion Green's function is: 



(Tc n (t)cl(Q)} 



(na - wi) 1+2 < 2S +!) (na + w<) 2 ( 2S + 1 ) 

^—ikpna 

{na + i;i) 1+2 ( 2 s+i) [n a _ w f) 2(2!+ ij 



(13) 



Considering now charge density wave (CDW), spin 
density wave (SDW), singlet superconductor (SS) and 
triplet superconductor (TS) correlations, it is easy to see 
that they all decay with the same exponent 2 + 2 s+i • 
These correlations are thus very strongly suppressed with 
respect to a free electron gas. Only for S — » 00, do we 
recover the exponents of the free electron gas for the con- 
ventional order parameters. If we specialize to N = 2, 
we obtain the same exponents as in Rcf . p5j . 

This leads us to investigate the presence of compos- 
ite order p9| , |l3| ]: the composite Charge Density Wave 

(c-CDW) order parameter O c -cdw — n(x).OsDw{x) 
and the composite singlet order parameter C _5 = 

n(x).OTs{x). 

It is easy to show that the composite Charge Density 
Wave order parameter can be expressed at the fixed point 



<J c -cdw ~ e 0i/2 



(14) 



where the field 0j/2° h a s conformal weight (j 
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4(25+2) ' 4 — 4(25+2) )• Similarly, for the composite sin- 
glet order, one has: 



Oc-s - e 



,1/2,0 
01/2 ■ 



(15) 



and both correlations decay with the same exponent, 

1 



{Ocomposite{n)0\ omposite {n )) 



7' Z 2S + 2 



(16) 



Therefore, for any S, composite order parameters are 
dominant. This can be understood in a simple way: elec- 
trons are tightly bound to local moments, so that com- 
posite correlations, that involve coherent motion of local 
moments excitations and electrons, decay more slowly 
than conventional ones. For S — > 00, conventional and 
composite order parameters become degenerate. This 
can be understood by remarking that in this limit, the 
local moments become classical and acquire a non-zero 
average so that composite and conventional order param- 
eters become identical. 

All our discussion up to now was concerned with the 
coupling of an integrable spin S chain with fermions. We 
now turn to the effect of the perturbation A$^$]/ that 
restores the behavior of the Heisenberg spin S chain |21| . 
If we assume that this term is small, i. e. that the Kondo 
coupling is larger than the largest gap induced by the 
perturbation for zero Kondo coupling, we can analyze its 
effect by determining whether it is relevant or irrelevant 
at the fixed point. The case of the spin 1 chain is special 
since Kac-Moody selection rules prevent the appearance 
of a primary operator of spin 1 in the fixed point theory. 
Generalizing (^) to include also non primary operators, 

a simple calculation shows that $^ — > Jr£' l , where Jr 
is the SU (2)i current, and £' L is a Majorana fermion of 
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the Ising (Mi minimal) model. This result can also be 
obtained from the Toulouse limit solution [^|[30|. As a 
result, the mass term becomes a term of dimension 3 at 
the fixed point and is irrelevant. We therefore expect that 
the underscreened regime we have found will be present if 
Xk exceeds the Haldane gap of the isolated spin 1 chain. 
For a spin S > 1, the operator <I>^ \ s becomes at the 



fixed point 



AO l $(!) 



This operator has dimension 

2S+1 ~ 2S+2 < 2 an d is therefore relevant at the fixed 
point. It is then likely that for S > 1 the models that we 
have described flow when perturbed to the trivial Kondo- 
Heisenberg fixed point even for a strong Kondo coupling. 

We may also examine what happens when two Kondo 
chains are coupled together to form a Kondo ladder. 
There are several ways to do it. One can couple 
the local moments of the Kondo chains by a Heisen- 
berg coupling Aj_ JY £>i,i-Si,2- At the fixed point, 
this gives rise to relevant terms of scaling dimension 
(25+if(25'+2) ' ^hat induce a spin gap. Second, one 
could consider an exchange interaction between the elec- 
trons, KeY / i c i,a,i^o ! ,0C i> i3 tl .cl, yt2 a r y > sCi t s ! 2. It is easily 
seen that such an interaction only leads to irrelevant 
terms, and does not affect our fixed point. Finally, 
one could consider interchain hopping of the electrons 

Si( c l o- i c i,o-.2 +h. c). At the fixed point, this leads to 
an operator of dimension 1 + 2 s+i ' re l evan t for S > 3/2. 
For 5=1, this term is marginal but it generates a rele- 
vant RKKY interaction between the local moments that 
destabilizes the chiral fixed point. 

We have obtained a physical picture of the one- 
dimensional underscreened Kondo Lattice. The forma- 
tion of a chiral non-Fermi Liquid results in strong an- 
tiferromagnetic fluctuations accompanied with compos- 
ite pairing. This picture is reminiscent of the situation 
that obtains in three dimensional heavy fermion systems 
||. It would be very worthwhile to try to determine if 
some analog of a chiral non-Fermi liquid can be found 
in higher dimensional Kondo Lattice models. A good 
starting point would be an array of Kondo- underscreened 
chains coupled by interchain hopping. 
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